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Abstract In this paper, we investigate various f{R)- 
brane models and compare their gravitational resonance 
structures with the corresponding general relativity (GR)- 
branes. Starting from some known GR-brane solutions, 
we derive thick /(i?)-brane solutions such that the met¬ 
ric, scalar field, and scalar potential coincide with those 
of the corresponding GR-branes. We find that for branes 
generated by a single or several canonical scalar fields, 
there is no obvious distinction between the GR-branes 
and corresponding /(i?)-branes in terms of gravitational 
resonance structure. Then we discuss the branes gener¬ 
ated by K-fields. In this case, there could exist huge 
differences between GR-branes and /(i?)-branes. 


1 lutroductiou 

Since general relativity (GR) was established by Ein¬ 
stein in 1915, a wide range of new theories of gravity 
have been proposed in the past 100 years. The sim¬ 
plest way to modify GR is to enlarge the dimensions of 
the space-time. The Kaluza-Klein (KK) theory as one 
of the pioneering five-dimensional space-time theories 
was proposed to unify electromagnetism with GR in 
the 1920s by Kaluza and Klein m- In order to gener¬ 
ate effective four-dimensional theories, the extra spatial 
dimensions in the KK types of theories are assumed to 
be compacted to the Planck scale. About 60 years later, 
physicists found that if our world is some kind of four¬ 
dimensional domain wall in a higher-dimensional space- 
time, then the matter fields can be dynamically trapped 
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on the wall. As a consequence, our world is effectively 
four-dimensional, even for infinite extra dimensions [3J 
mOiniE]- String theory also predicts the existence of 
TeV scale extra dimensions [5]. 

But what really triggered the modern revolution of 
extra dimension theories is the work done by Arkani- 
Hamed, Dimopoulos, and Dvali (ADD) [9], which pro¬ 
vides an alternative solution to the gauge hierarchy 
problem. The idea of ADD model is that all the matter 
fields are confined on the four-dimensional brane em¬ 
bedded in a higher-dimensional space-time, while grav¬ 
ity can propagate in the bulk. In ADD model, the fun¬ 
damental scale is assumed to be the electroweak scale 
M* = mEw ~ 1 TeV, and the large hierarchy between 
the Planck scale Mp\ and the fundamental scale is as¬ 
cribed to the large volume of the extra dimensions. 
Due to the compactification of the extra dimensions, 
the four-dimensional Newton law appears at a distance 
much larger than the radius of the extra dimensions. 

Soon after the ADD model, Randall and Sundrum 
(RS) found that by using a nonfactorizable warped ge¬ 
ometry, one can also solve the hierarchy problem [lOj . 
More interestingly, they have shown that it is possi¬ 
ble to generate the four-dimensional Newtonian law at 
large distances in a warped space-time even when the 
extra dimension is infinitely large m- Now, these two 
models are known as the RS-I and the RS-II models, 
respectively. In the RS-II model, the four-dimensional 
Newtonian law is generated by a gravitational zero mode 
localized near the brane. Later, Gregory, Rubakov, and 
Sibiryakov (the GRS model) [12] have shown that four¬ 
dimensional gravity can also be reproduced at an inter¬ 
mediate scale even without a localizable gravitational 
zero mode. In fact, four-dimensional Newtonian gravity 
is generated by a quasilocalized graviton (or gravita- 
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tional resonance) in the GRS model [HlfH] . The GRS 
model also shed light on solving the cosmological con¬ 
stant problem But on the other hand, it suffers 

from problems such as violating the positivity of energy 
and containing ghost modes [laiiaiiiiiiTiiiaiis]. 

Another interesting and widely studied model which 
reproduces four-dimensional Newton’s law at interme- 


been applied in a wide range of cosmology studies m 
(see for reviews, and [531 

for f{R) braneworld). 
The graviton equations for a class of thick /(i?)-brane 
models were derived in Ref. m- It is shown that the 
graviton equation is determined by both the background 
solution and the function f{R). Using the results of 


diate scale is the DGP model proposed by Dvali, GabadadzeRef. [77], some of us find a series of gravitational reso¬ 


und Porrati [50]. The DGP model offers us a new way 
to understand the acceleration of the universe ED- It 
is interesting to note that in higher-dimensional real¬ 
ization of the DGP model, the four-dimensional New¬ 
tonian gravity is reproduced by massive graviton res¬ 
onances [22l[23l[24] . In addition to the GRS and DGP 
models where the bulk is asymptotically flat, some works 
have considered graviton resonances in asymptotical 
AdSs space Gravitational resonance also 

appears in other topic of gravity research. For instance, 
in black hole perturbation theory it is possible to find 
gravitational resonances which are known as the quasi¬ 
normal modes m- No doubt, gravitational resonances 
deserve further consideration. 

In this paper we focus on the appearance of grav¬ 
itational resonances in thick RS-II type brane models, 
where the brane can be generated by, for example, a 
single canonical scalar field in an AdSs space [sniiniisa 
[55] (more complicated thick brane models can be found 
in ED). The possibility of the existence of gravitational 
resonances in thick brane models was first pointed out 
in Ref. [31] . But unfortunately, the solution therein does 
not support any gravitational resonance. In thick brane 
models constructed in general relativity, gravitational 
resonances can be found in cases with either a single 
canonical background scalar field [311111] , or two canoni¬ 
cal scalar fields (the Bloch brane [37]) [551115] . or a single 
noncanonical scalar field (also known as K-field [iOlHTl 
132] ; see [i51imi551l551IT7llT5] for thick branes generated 
by K-fields) [55] . 

However, the spectrum of graviton is determined 
by the equation of the gravitational mode (the spin- 
2 transverse and traceless tensor mode of the metric 
perturbations my In general relativity, the equation 
of the gravitational mode is independent from the back¬ 
ground scalar fields and only depends on the warp fac¬ 
tor. But this is not always true in modified gravity the¬ 
ories. Thus, in order to find gravitational resonances, 
it is interesting to go beyond general relativity. For ex¬ 
ample, it is possible to find gravitational resonances in 
thick brane models in scalar-tensor gravity [50] . 

Another simply modified gravity is f{R) gravity, 
where the Lagrangian of gravity is an arbitrary func¬ 
tion of the scalar curvature. f{R) gravity was first in¬ 
troduced in cosmology in 1970 and nowadays it has 


nances in Ref. [75], where the gravity Lagrangian den¬ 
sity is /(i?) = R + aR^, and the brane is generated by a 
single canonical scalar field. However, there exist some 
singularities in the model of [75] . 

The aim of this paper is to construct thick f{R)- 
brane models which support gravitational resonances 
but are free of singularities. By starting with some well- 
studied GR-brane solutions, we construct new f{R)- 
brane solutions (either analytically or numerically) such 
that the metric and the scalar configurations remain 
the same. Then we compare the gravitational resonance 
spectra between our /(i?)-branes and their correspond¬ 
ing GR-branes. Since gravitons satisfy different equa¬ 
tions in different gravity theories, there might be some 
differences in the resonance spectra. 

Note that in addition to the tensor mode, there 
are also spin-1 vector modes and spin-0 scalar modes. 
These modes evolve independently, no matter the thick 
brane is constructed in general relativity [751155115711551 
[531, or in /(R) gravity [5H1I55] . The gravitational reso¬ 
nances only relate to the tensor sector. Usually, the vec¬ 
tor modes only contain an nonnormalizable zero mode 
and therefore can be neglected when one only consid¬ 
ers static sources [751l801l831l84] . The scalar sector will 
also contribute to the four-dimensional newtonian grav¬ 
ity [55]. However, for /(i?)-brane generated by a back¬ 
ground scalar field, the scalar sector of the linear per¬ 
turbations has two degrees of freedom after fixing the 
gauge. In general, these two scalar modes are governed 
by two coupled second-order differential equations, which 
are too complicated to be considered here and deserve 
an independent work elsewhere. 

In our first /(i?)-brane model, we briefly discuss the 
similar model considered in Ref. [75] and add the con¬ 
dition fa = df{R)/dR > 0, which guarantees that the 
graviton is not a ghost and the normalization of the 
zero mode. In this model, we find that there is nothing 
special about gravitational resonances in /(i?)-brane 
model in contrast to GR-brane. In order to get a differ¬ 
ent behavior, naturally, one can generalize this model to 
the case with multiple scalar fields. Based on this mo¬ 
tivation, we study the brane generated by two scalar 
fields in second model and we obtain similar results as 
the case with a single scalar field. Our analysis shows 
that for the brane generated by canonical scalar fields 
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(no matter how many they are), there is no striking dif¬ 
ference between /(i?)-brane and GR-brane in terms of 
gravitational resonance. Under the circumstances, non- 
canonical scalar fields may be a breakthrough. In our 
third /(i?)-brane model, indeed, we find that even for 
the case with only one single noncanonical scalar field, 
gravitational resonances on /(i?)-brane can be signifi¬ 
cantly different from the GR-brane. 

Therefore our work is organized as follows. In Sect. [21 
we review the tensor perturbations of /(i?)-brane and 
the condition of normalizable zero mode. In sect. |3l we 
construct some thick /(i?)-brane models in cases with, 
first, a single canonical scalar, and then two canonical 
scalars, and finally a single noncanonical scalar field or 
K-field. These /(i?)-brane models are all based on the 
corresponding GR-models which have been extensively 
studied in the literature. Then we compare, via both 
analytical and numerical methods, gravitational reso¬ 
nances between our /(i?)-branes and the correspond¬ 
ing GR-branes. In Sect. 01 we analyze how the f{R)- 
brane models impact gravitational resonances. Finally, 
in Sect. 01 we come to the conclusions and discussions. 


2 Tensor perturbations of /(ii)-brane 

Now we consider the five-dimensional /(R)-brane model. 
Our action is taken in the general form 

S = I d^xdy^ , (1) 

where f{R) is an arbitrary function of the scalar curva¬ 
ture R and = SttGs with G 5 the five-dimensional 
Newton constant. For convenience, we take K 5 = 2. 
Capital Latin indices M, N,... = 0,1, 2, 3,4 denote the 
bulk coordinates, and Greek indices iz,... = 0,1, 2, 3 
denote the brane coordinates. The coordinate y = 
stands for the extra dimension. (j)i is the ith scalar field, 
and its kinetic term Xi is given by 

Xi =(j)idM(t)i- ( 2 ) 

We consider a static flat brane, for which the line 
element can be written as 

ds^ = a^{y)rii_,,,dx^dx'' + dy^, (3) 

where a^{y) = is the warp factor, and 77 ^;, is 

the induced metric on the brane. For this background 
space-time, the scalar fields only depend on the extra 
dimension, i.e., (j)i = fpiiy)- The scalar curvature is given 

by 

i?(y) =-4(5A'2-P2yl"). (4) 


In this paper the prime denotes the derivative with re¬ 
spect to y. The Einstein equations and the equations of 
motion for the scalar helds are 

f + 2fR (4^2 + yl") - - 2/"+4L = 0, (5a) 

f+8fR{A"+A'^)-8rRA'+dJ2^xA?+^L = 0, (5b) 

i 

and 


4'i{Lxi +2XiLxiXi) + L^,^=2XiLxi<i)i—dLxiCt^'iA', ( 6 ) 

respectively, where Lxi = dLjdXi, L^. = dL/d(j)i, and 
fn = df{R)/dR. 

In this paper, we only consider the minimal coupling 
between (pi and Xi. Therefore, the Lagrangian density 
L((/)i, W) is given by 

L{(P„Xi) = Y,[F{X,)-V{cP,)], (7) 

i 

where F{Xi) is an arbitrary function of Xi, and V{(j)i) 
is the scalar potential of the scalar field (pi. 

Next, we consider the tensor perturbations of the 
background metric ([3|): 

ds^ = + h^^)dx^dx'' + dy^, ( 8 ) 


where satishes the transverse and traceless (TT) 
conditions = 0 = d^t.h'p. With the coordinate 

transformation dz = a~^dy, the perturbed Einstein 
equations read ua 


+ + 

V a JR ^ 


□ (4) 


hu]., — 0 , 


(9) 


where Note that although in Ref. [77] 

only one scalar field with canonical dynamics was con¬ 
sidered, Eq. remains valid for our case because the 
TT tensor modes are decoupled from other modes of 
the perturbations. Then, with the KK decomposition 

m 

hfj,^{xP,z) = efj,^{xP)a~^/'^fP^^^'^tp{z), ( 10 ) 

where e^vix^) satisfies the TT conditions = 0 = 

Eq. dH) can be simplified as a Schrddinger-like 
equation for 7 />(z) [77] : 


[-dl + W{z)] ip{z) = m^iP{z). 


( 11 ) 


Here the effective potential W{z) is given by [77] 

3 {d^aY 3^ 8d:,ad^fR 

A a2 2 a 2 ' 

^ 1 dlfR 
2 


W{z)=\ 


afn 


1 

fR ) 


Jr 


( 12 ) 


In the next section we will discuss the effective potential 
W(z(y)) in the coordinate y. Using the relation dz = 
ady, Eq. (02]) can be rewritten in the coordinate y: 

W{z{y)) = -k ^(a'- 


4V fR ) 


„ 2>aa' fL 

Taa")T—^ 

laa%+^ 

2 fR 


(13) 
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From the Schrodinger-like equation (ED, it is easy 
to give the solution of the gravitational zero mode 

= (14) 

where Nq is the normalization constant. To get four¬ 
dimensional Newtonian gravity on the brane, the zero 
mode should satisfy the following normalization condi¬ 
tion: 

/ CX) 

\^^°'>{z)\'^dz <oo. (15) 

-OO 

Apparently, in order to ensure the normalization of the 
zero mode, fn should satisfy Jr > 0 in the whole 
space, which also guarantees that the graviton is not 
a ghost E] • 


3 Gravity resonances in various /(i2)-brane 
models 

In this section we study braneworld solutions and grav¬ 
itational resonances in three types of thick /(i?)-brane 
models. In our first model, the brane is generated by a 
single canonical scalar field, which is one of the most 
widely studied thick brane models. Next, we consider 
the brane generated by two canonical scalar fields (known 
as the Bloch brane) in the second model. The main 
property of the Bloch brane is that the brane is split 
into two sub-branes. Finally, we study the brane gen¬ 
erated by a K-field, whose most important feature is 
that its kinetic term is noncanonical, which will pro¬ 
duce some different results compared to the previous 
two cases. 


3.1 /(i?)-brane model with L = X — V{<f)) 

First of all, we study the /(i?)-brane generated by a 
single canonical scalar held, whose Lagrangian density 
reads L = X — V{4>) = This kind 

of /(i?)-brane model is one of the simplest and the most 
widely studied models. Therefore, one reason why we 
study it is to review and illustrate the research pro¬ 
cesses of gravitational resonances on /(i?)-brane. The 
second and also more important reason is to give a di¬ 
rect comparison with our third model generated by a 
single noncanonical scalar held. The dynamical equa¬ 
tions m and dHD become 

f + 2fR (4A'2+A")-6/)^A'-2/"=2(<(>'2 + 2y), (16a) 

-8fR {A” + A'^)+8fRA' - f=2{^'^-2V), (16b) 

<(>"+4AV'=y^. (16c) 


For the case of f{R) = -R, we consider the solution 
given in m- 


(j){y) = arctan 


tanh 


A{y) = —6 In [cosh(fcy)], 


V{4’) = 




(1 — 4&) -I- (1 -I- 4&) cos 



(17a) 

(17b) 

(17c) 


where h and k are positive parameters. For convenience, 
we set y = ky in the next content, which is a dimen¬ 
sionless parameter. 

Now we construct the /(i?)-brane model that shares 
the same background solution (I17a|) of the GR-brane. 
Keeping the above background solution of the GR-brane 
unchanged, and incorporating Eqs. (I16al) and (I16bl) , the 
function Jr should satisfy 


- fR + fRA-8fRA = 2<p^>0, (18) 

where A = —b tanh(y), A = —bsech^{y), = |6sech^(y), 
and the dot denotes the derivative with respect to y. 

It is interesting to give the explicit expression for 
fR and f{R). The procedures are described as follows. 
First, we solve Eq. (TT^ to give the analytic solution for 
Jr as a function of y: 


//j=l-hasech2 {y) (tanh(y))-/3Q|._ (tanh(y)) 


(19) 


where a is an arbitrary constant, P and Q are, respec¬ 
tively, associated Legendre functions of the first and 
second kinds, /3 = Pr^{0)/Q%'^{0) is given by the con¬ 
dition /fl|y=o = 0, and K± = ^^/{b— 14)6 -|-1 ± 1/2. 
Second, we define R=^=Ab [(56 -I- 2)sech^(y) — 56], 
which is obtained from Eqs. (0]) and (I17bl) . to give the 
expression of y in terms of the curvature scalar R: 


y = arccosh 


f 2V26-t562 ' 

\VmA~+R, 


( 20 ) 


Thirdly, defining f=f/k'^, which results in = /r, 
and substituting y{R) into (fTO|) yield the expression of 
as a function of R. Einally, integrating fj^ with re¬ 
spect to R gives the expression of f{R). In the following, 
we will first consider some special values of the param¬ 
eter 6, which will result in simple solution of /(R), and 
then general 6. 


3.1.1 6= 1 

For the case of 6 = 1, we obtain the following solution: 
f = 1 + ai cos (n{R )), (21) 
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Ci 

Value 

Cl 

240936850868416 

149442323051763 

C2 

38202658243 

16604702561307 

C3 

2065663433 

4605037510335808 

C4 

18231247 

3831391208599392256 

C5 

351657 

3187717485554694356992 

C6 

6507 

2652180947981505705017344 

C7 

459 

15446301841044289226021011456 


Table 1 The values of the parameters Ci in Eq. 


where 'H{R) is given by 


( 22 ) 


Note that satisfies /^|y=o = 1 + cti. Then the func¬ 
tion f[R) can be integrated from Eq. (1^ as 


yv 2o'i r 

f{R)=R + -y [2(i? -h 6) cos 


( r { r )) 

Y^3(.R-8)(i?-h20) sin [u^R^ 


(23) 


where the integration constant is absorbed into the scalar 
potential E (<()). 


3.1.2 6 = f 

For the special case of 6 = ^, we have 

7 

f\R)=R + a2j2(^nR^’ (24) 

n—1 

where the values of the coefficients Ci are given in Ta¬ 
ble [TJ The advantage of taking those values is that we 
can guarantee /^|y=o = ^ + ct 2 , which is similar to the 
case 6=1. 


3.1.3 The general b 


For general 6, Eqs. m and (EHl) give the solution of 
/^, which is a complex expression. In this case, we can 
give the analytic f{R) from Eqs. (IT^ . (1^ . and (llbaj) 
or (flBbll . which is given by 


f{R)=R+a 




2-h56 

-Aib‘^-2bK+)T'\p^^{W) 


+pw (QT.m - Q^lm) 



(25) 


where the parameters a, j3^ and K± are those appearing 
in m and 0 and P are given by 


p= ^ 1 - 6 ) 2 , 

^ _ \/2062 + R 

~ 2V26 -H 562 ■ 


(26) 

(27) 


To investigate the resonant modes of gravity we re¬ 
call the definition of the relative probability [5S] : 


P{m^) 




(28) 


where ’4>{z) is the solution of Eq. (fTTll . 2zb is approxi¬ 
mately the width of the brane, and Zmax = 10z{,. Here 
|'i/'(2)P can be explained as the probability density [SSJ 
I86j . A resonant mode with mass exists, if the rela¬ 
tive probability P{m3) has a peak around m = m„ and 
this peak has a full width at half maximum (FWHM). 
So, the number of the peaks that have FWHM corre¬ 
sponds to the number of the resonant modes. For the 
sake of simplicity, we impose the following conditions 
on tp{z)-. 


'0even(O)=l, l9z'0even(O) = 0; (29a) 

^odd(0)=0, 9,V^odd(0) = 1, (29b) 


where tpeven and '(/’odd correspond to the even and odd 
parity modes of respectively. Since there is no es¬ 
sential difference between the above conditions for the 
relative probability P{nR), we will only display our re¬ 
sults with one of the conditions (12911 in the following 
examples. For the sake of brevity, we discuss gravita¬ 
tional resonance in the coordinate y, but set y = y, 
R = R and / = /, which is equivalent to setting k = 1 
in the coordinate y. 

Note that in this /(i?)-brane model there are two 
primary parameters 6 and a. The parameter 6 deter¬ 
mines the brane solution and a affects the range of 
values of the function f{R). Therefore, they may char¬ 
acterize the gravitational resonances that will be dis¬ 
cussed in the following. It can be seen that the expres¬ 
sion of f{R) in ((Ml) for general 6 is very complex while 
it becomes simple for two special values: 6 = 1 and 
6 = 52/3. We find that the shapes of f{R{y)) as the 
functions of y for the two special and general values of 
6 are similar. Therefore, in the following calculations, 
we only need consider the two special values, for which 
the expressions of f{R) and the effective potential W(y) 
are simple. As for the parameter a, it characterizes the 
deviation of the /(i?)-brane from the GR-brane. This 
can be seen from (UK). 

For the case of 6 = 1, we substitute (l23ll into o 
to solve Schrodinger-like Eq m with three typical val¬ 
ues of ai: ai =0 (GR-brane), oi = 100 (/(i?)-brane 
with fn > 1), and ai = —0.13 (/(i?)-brane with 0 < 
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(b) ai = 0, -0.13 
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(c) ai =0,100,-0.13 


(d) ai = 0, 100, -0.13 


Fig. 1 Plots of the function fji, effective potential W{z{y)), and relative probability P{rn?) of the odd KK modes for the 
single scalar brane with standard kinetic term with 6 = 1. In the figures, the GR-brane (ai = 0) is denoted by the red 
solid lines, the /(i?)-branes with ai = 100 and ct\ = —0.13 are donated by the black dotted lines and blue dot-dashed lines, 
respectively. (From now on, red solid lines always correspond to the GR-brane, black dotted lines correspond to fn > 1, blue 
dot-dashed lines correspond to fn < 1, and the rest of figures are denoted by black solid lines in all our plots). 


fii < 1) (see Fig. 1(a) 1(b) I. Note that the functions 
fniy) for cti = 100 and ai = —0.13 (both guarantee 
fniy) > 0 in the whole space) have already deviated 
from fniy) = 1 for ai = 0 (the GR-brane) as far as 
possible. For a > 100, the potential W{z{y)) can hardly 
change more compared with ai = 100. For ai < —0.13, 
fii may be negative at some points. The effective po¬ 
tentials W{z{y)) are shown in Fig. 1(c) Our calculation 
shows that for all the chosen values of ai (see Fig. |l(d)] ), 
there is no gravitational resonance. 


For the case oi b = the corresponding solution 
of f{R) is given by Eq. (1^ . Compared with the case of 
6=1, although the maximum of the potential W{z{y)) 
is enlarged, but no resonant KK mode is found. The cor¬ 
responding results are shown in Fig. [2] For other values 
of 6, the results are similar: the /(i?)-brane can hardly 
deviate from the corresponding GR-brane in terms of 
gravitational resonance and we do not find any gravi¬ 
tational resonance. 


3.2 /(i?)-Bloch brane model with 
L = Xi + X 2 — V (^ 1 , (j)2) 


Next let us discuss the Bloch brane, which is a kind 
of important brane with inner structure generated by 
two scalar helds. It was hrst obtained in Ref. m and 
then developed and discussed in Refs. [87113^1881I891IW)1 
. In this subsection, we explore gravitational res¬ 
onances on the Bloch brane in the f{R) gravity (called 
the /(i?)-Bloch brane) for several reasons. First, there 
exist gravitational resonances on the Bloch brane [33] . 
Therefore, we can compare the difference between the 
GR-Bloch brane and /(i?)-Bloch brane in the presence 
of gravitational resonances, which is expected to be a 
complement to the former /(i?)-brane model. Second, 
the Bloch brane is generated by two canonical scalar 
fields. Investigation of the Bloch brane in this subsec¬ 
tion could pave the way for Sect. 01 where we will gener¬ 
alize the conclusion of the first two models to the cases 
with multiple canonical scalar fields. 
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fRiy) 
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Fig. 2 Plots of the function fR{y), effective potential W{z{y)), and relative probability P{m?) of the odd KK modes for the 
canonical single scalar brane with 6 = 52/3 and 012 = 0, 100, —0.012, respectively. 


The Lagrangian density © for the Bloch brane is 
given hy L = -\d^4>dM(t> - ^d^x^MX - V{4),x), 
where (j) Eind x two interacted real scalar fields. 

The dynamical equations ([S]) and © are modified 
to 


/ + 2/fi (44'2+4") - x'2^ 

-/ - 8/« {A" + A'^) + x'"-2y),(30b) 

(j)"+AA^<t)'= V^, (30c) 

X" + 44Y = K,. (30d) 


For the case of & > 2d > 0, the general symmetric 
Bloch-brane solution reads [371157] 


4>{y)=v tanh(2dtit/), 
Xiy)='v\h—^ sech(2dvy), 


A(y )=—- (b— 3d) tanh^ (‘^dvy) — 2b In cosh(2c?t;y) 
9a L 


(32a) 

(32b) 

(32c) 


Analogously, selecting certain values of the param¬ 
eters V and b, we can calculate the explicit expression 
of f{R). For V = \f 3 J 2 and b = 3d, we have 


For f{R) = R, there are many symmetric and asym¬ 
metric solutions I371I571I57] . We consider one of the so¬ 
lutions given in Refs. [371157], where the scalar potential 
V{(p) takes the following form: 


V[P) = 


(iv^-lp^-dx^) +4dYV 
1 . ^ 2 


- o - dpx^ 


(31) 


/k = 1 + 7 C 0 S iy{R)), 


(33) 


where 7 is a constant, R = 24 ^^ [7 sech^ (-x/Ody) — 5], 
and y{R) reads 


y{R) = V31n 


Vi?-48d2-p 


2V42d 

Therefore, integrating the function fn over R we get 
27 


(34) 


f{R)=R + y [\/3(i? - 48d2)(i20d2 -p R) sAvy{R) 
+2{36d‘^ + R)cosy{R)]. (35) 
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Fig. 3 Plots of the function fR{y), effective potential W{z{y)), and relative probability P{m^) of odd resonance KK modes 
on the Bloch brane with d = 1, = 1, cq = —2 — 7.0 x 10“^®, and fjily—Q = 1, 50, 0.92, respectively. 


Note that when d = y this solution is identical with 
Eq. (ESI). 


For the case oi b = d, the authors of Ref. [57] ob¬ 
tained the degenerate Bloch-brane solution: 


(t>{y)=t 


x{y)= 


Vc 




\/co- 




.v4 


g — 4 sinh(2dz;j/) 

- 4cosh(2dtiy) — cq 

2v 

4cosh(2dtiy) — cq 

a/c§ - 4 - Cq 

- 4cosh(2duj/) — Co 


4v^/9 


(36a) 

(36b) 


xexp 


xexp 


4t;2 eg - 4 - coi/cg - 4^ 

9 (\/cg - 4 - Co) 

4t;^ ^cg — 4—Co a/ eg — 4 cosh(2dvy)^ 
9 ^i/cg —4cosh(2dut/) —co^ 


(36c) 


where cq < —2 . 
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Next, we take the above Bloch-brane solution (ESI) 
as our background solution of /(i?)-brane. Using the 
similar process used in the previous section, the func¬ 
tion fji should satisfy 

- /fl + f'R^' - ifR^" = + x'") > 0. (37) 


Since it is difficult to give the analytical solution of Jr, 
we can obtain Jr (also f{R)) and the potential W{z{y)) 
by means of numerical methods. Because f{R{y)) (and 
hence Jr) is an even function of the extra dimension 
y, we have /)j|y=o = 0. To construct a numerical Jr 
satisfying Eq. Eli), we need one more initial condition, 
i.e., the value of /ij|y=o- With the adjustable parameter 
/fl|y= 0 j we can construct a different numerical Jr. 

Here, we consider the following typical set of param¬ 
eters reducing a double-brane: a = (i=l, v = 2, and 
Co = —2C7.0 X 10“^®, which was used in Ref. EH] for 
the GR-brane model. It is convenient to compare the 
results of the /(i?)-brane and the GR-brane. 

The effective potential W{z{y)) and the relative prob¬ 
ability P{m?) of the gravitational KK modes are shown, 
respectively, in Fig. 3(c) 3(d)[ We find that there are 


five resonance peaks in Fig. |3(d) for the GR-brane. 
For the /(i?)-brane model, the potential W{z(y)) (see 
Fig. |3(c)[ ) almost reaches the maximum deviation from 
the GR-brane when = 50 and 0.92 (the nu¬ 

merical solutions of /r are shown in Fig. 3(a 3 O. 
The corresponding resonance mass spectra are shown 
in Fig. |3(d)[ from which we can see that there exist four 
peaks for /i^|y=o = 50 and five peaks for /fl|y=o = 0.92, 
respectively. Therefore, there is no obvious difference 
between these spectra. For other values of the parame¬ 
ters we also have the same conclusion. 


3.3 /(i?)-brane model with L = X — XX'^ — V{(j)) 


dynamical equations © and © are reduced to 
f+2fR (4H'2+W')-6/;jH'- 2/"=2 (^0'2+iA</>'4+2U^ , 

-/- 8/fl (H"+H'2) + 8/;jH'=2 ^ (38) 


For f{R) = i?, the authors in Ref. [57] obtained one 
analytic solution. For general /(R), in the vast majority 
of cases, the dynamical equations remain hard to solve 
except in some specific cases (for example, if f{R) is a 
polynomial). As an example, we consider the following 
scalar field and warp factor: 


tpiy) = w arctan 



A(y) = -ln [cosh( 2 /)]. 


(39) 

(40) 


For the case of w = 1, we get the solution 


f{R) = R + 


A 


3136 




(29A - 980) 
n76 




(41) 


where R = 28 sech^{y) — 20 and A > — ^ to guaran¬ 
tee Jr > 0 . 

If we require the coefficient of i? to be 1, then we 
obtain 


, /±14Vl74A-b49-98 r , 

(j){y)=\l --arctan j^tanh 

A{y)=- In [cosh(y)], 
fm R (87AT7V174A + 49 + 49) - 
' 6728A 


(I) 


R\ 


(42) 

(43) 

(44) 


For -I- in Eq. (l42|) . we require — ^ < A < | and — 
corresponds to — < A < — ^. This solution was also 

obtained in Ref. |76j recently with another method. 

Here, to study the solution of (j) with more abun¬ 
dant structure, we reconstruct one analytic solution of 
(f) given as 


In the previous two subsections, the /(i?)-branes are 
all generated by one or two canonical scalar fields. For 
comparison, we try to investigate /(i?)-brane generated 
by noncanonical scalar fields (K-fields) in this subsec¬ 
tion. 

K-hleds were introduced to study inflation for the 
first time gnuniiia and then were extended to many 
fields. In brane world models, branes generated by K- 
fields is generally known as K-branes, which might present 
some new properties (such as the localization of bulk 
fermions and gravitons m) compared with the corre¬ 
sponding standard branes. Therefore, we expect there 
are some new results of gravitational resonances on K- 
branes. 

We consider the simple case of one K-field with the 
Lagrangian density L = X — AA^ — V{(j)) (A 7 ^ 0). The 


(j) = barcsin [tanh(?/)] -|- c?tanh(?/)sech(y), (45) 

where b and d are free parameters. Therefore, the warp 
factor is determined by 

2X(I)"^ + 3A" -t 20'2 ^ 0 , (46) 

with the boundary conditions A( 0 ) = 0 and A'( 0 ) = 0 . 
Meanwhile, /r should satisfy 

2A<(.'4 + 3/flA" + /" - f^A^ + 2(/.'2 = 0 . (47) 


The value of fR\y=Q can be arbitrary so long as the 
corresponding Jr satisfies Eq. (|57|) and guarantees Jr > 
0 in the whole space. 

If the parameter A is positive, the solution of the 
wrap factor will be lumplike and the constraint condi¬ 
tion (1571) will be the same as Eqs. (fl^ and (1571) . Our re¬ 
sults show that it has no breakthrough comparing with 
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0(y) 


0(y) 





Fig. 4 Plots of the scalar field (j> given by Eq. I(i5l). The figures are, from left to right, kink-like solution, non-monotonic kink 
solution, and double kink-like solution. 


W(y) P 




Fig. 5 Plots of the effective potential W{z(y)) and relative probability P(m?) of the odd resonance KK modes on the brane 
generated by a K-field with 6 = 2, d = 2, A = 0.6, and fji\y—Q = 1,100. 


the former cases discussed in the previous two subsec¬ 
tions. 

For A > 0, we consider three types of solutions of the 
scalar field (see Fig. 0]). Because when 0 < < 1 

the range for //j|y=o is so narrow, we only consider the 
case /ij|y=o > 1 in this /(i?)-brane model. We find that 
when the scalar field (j> is taken as the form of single kink 
or non-monotonic kink, the results are the same as the 
brane generated by a canonical scalar field (see Figs.[5l 
©. For the case of double kink, the result is similar to 
the Bloch brane (see Fig. 0. 

For A < 0, we also set three similar groups of param¬ 
eters (see Fig. and all of them have lumplike solutions 
for warp factor. Here we also only consider the values 
of parameters which result in /_R|y=o ^ 1 in /(i?)-brane 
model. Obviously, the consequences (see Figs. Enniiiii) 
are analogous to the case A > 0. 

It seems that we can draw the conclusion that the 
difference of the gravitational resonances is small be¬ 
tween the /(i?)-brane and GR-brane when the warp 
factor is monotonic in the area y £ [0, -boo). 


Therefore, we consider another interesting scenario, 
i.e., the warp factor is a non-monotonic fnnction of y 
in the area y £ [0,-|-oo). The parameters are taken as 
6 = 2, d = —2, and A = —0.6 (see Fig. [121). The effective 
potential W{z{y)) and relative probability P{m?) of the 
GR-brane are shown in Fig. 13(a) 13(c) Apparently, 
there is no gravitational resonance on the GR-brane. 
Setting the same parameters for the /(i?)-brane model, 
it can be seen that when fii\y=o = 0.11, in compari¬ 
son to the GR-brane the effective potential and relative 
probability change greatly (see Fig. 13(b)[ 13(c)). 

For the set of parameters 6 = 2, d = —2, A = —0.6 , 
and f}i\y^Q = 0.106, we can see that fn almost vanishes 
at 2 / = ±1 (see Fig. 14(a) I. We display the effective 
potential W{z(y)) in Fig. 14(b) and relative probability 
P{m?) in Fig. 14(c) to prepare for the next section. 


4 The impact of /(i?)-brane 

Now, we analyze how /(i?)-brane model impacts gravi¬ 
tational resonances. Focnsing on Eq. m, the different 
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W(y) 




m 


(a) fR\y=o = 1,100 


(b) /H|y=o = 1,100 


Fig. 6 Plots of the effective potential W{z{y)) and relative probability of odd resonance KK modes on the brane 

generated by a K-field with b = l,d = —2, A = 0.6, and fR\y—Q = 1, 100. 

W(y) P 




m 


(a) fR\y=o = 1,100 


(b) /H|y=o = 1,100 


Fig. 7 Plots of the effective potential W{z(y)) and relative probability P(m?) of the odd resonance KK modes on the brane 
generated by K-field with b = 6,d = —6, A = 0.6, and fR\y—Q = 1, 100. 


parts of the potential W{z{y)) between the GR-brane 
and /(i?)-brane are denoted by AW: 

fk UafnY I aa'fn + . 43 ) 

2 Jr Jr J 2 /r 'VI 

First of all, we analyze the model with only one 
canonical scalar field. In this case, we have for the GR- 
brane 

= (49) 

which results in a lumplike solution of the wrap factor. 

For /(i?)-brane, we have Eq. (fT51) . Fixing the so¬ 
lution of the warp factor, if we expect some huge dif¬ 
ferences of the potentials W{z{y)) between /(i?)-brane 
and GR-brane, the value of jf must be large 

enough at some points to form deep wells or high bar¬ 
riers. Gonsidering that f{R) should be continuous and 
satisfy fR>0, the most promising approach may be 
that /r approaches zero at some locations. 

Because both f{R) and /r are even functions of 
the extra dimension y and, generally, f{R) should be a 


constant at infinity, we need only consider the situation 
that Jr approaches a very small value at one or two 
locations (more locations makes no difference). 

From Eq. dm) we know that AW is also an even 
function. We assume AW has two infinitely deep po¬ 
tential wells (the same analysis is appropriate for bar¬ 
riers) symmetrically situated at y = iyo- For simplify, 
we consider the case that these points are minima of 
fR{y). So, we have = 0 and /)( > 0 at these points. 
We plug /){ = 0 into Eq. dm at y = ±yo, then we have 

-3kl"/K > f'A- 

Obviously, the first two terms of AW in Eq. (HRl) 
vanish at y = ±yo. Since —iA"J r > fA, the ranges 
of AW at y = ±yo are both 0 < AW = < 

which is finite. So, it is hard for the 
potential W{z{y)) to have deep potential wells or bar¬ 
riers around y = ±yo even if fji\y^±y^ approach zero. 
The same analysis also applies to the situation that 
W{z{y)) is supposed to have one infinite potential well 
at y = 0. 
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0(y) 0(y) 0(y) 





Fig. 8 Plots of the scalar field (j> given by Eq. I(i5l). The figures are, from left to right, a kink-like solution, a non-monotonic 
kink solution, and a double kink-like solution. 

W(y) P 




(a) /h|j,=o = 1.50 


(b) /h|i,=o = 1,50 


Fig. 9 Plots of the effective potential W{z(y)) and relative probability P(m?) of the odd resonance KK modes on the brane 
generated by K-field with 6 = 1, d = 0.3, A = —0.01, and fji\y—Q = 1, 50. 


As for the Bloch brane, although the solutions of 
the warp factor and scalar fields are not the same as 
the previous model, the restrictions on the warp factor 
and f{R) are similar (see Eq. (I^ h So, we also have 
the same conclusion. 


Finally, we discuss why the third /(i?)-brane model 
has a strong influence on the potential W{z{y)). Taking 
f[R) = i? in Eq. (H71) . we find that although A is a small 
parameter, it can lead to a non-monotonic warp factor 
from y = Q to y = -boo. Therefore, the shape of the 
potential W{z{y)) of the /(i?)-brane becomes compli¬ 
cated. Furthermore, obviously when A < 0, it is possible 
that the function fn^y) approaches zero at some points. 
In Fig. 14(a) we see that when fniy) gets close to zero, 
there are two high barriers at the corresponding points 
(Fig. 14(b) I but the GR-brane never has such a kind of 
feature. 


5 Conclusion and discussion 

In this work we discussed gravitational resonances in 
various /(i?)-brane models. For the branes generated 
by one or two canonical scalar fields, we found that the 
/(i?)-brane and GR-brane are similar from the view 
point of gravitational resonance. Then we investigated 
the brane generated by a single K-field. For a positive 
A (or a negative A with the warp factor varying mono- 
tonically), the structure of gravitational resonance of 
/(i?)-brane does not deviate much from the GR-brane. 
This result is similar to the case with the standard ki¬ 
netic term. When the parameter A is negative and the 
wrap factor is non-monotonic, the structure of the grav¬ 
itational resonance of the /(i?)-brane can seriously de¬ 
viate from the GR-brane. The degree of deviation relies 
mostly on the value of the parameter fnly^o- 

Finally, we analyzed why the /(i?)-brane and GR- 
brane have similar gravitational resonance structures 
when the background scalar field is canonical. Our con¬ 
clusion is that the similarity is due to the limitations 
fR>0 and Eq. ([THl), which constrain strictly the po- 
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W(y) 


P 




(a) /h|j,=o = 1.50 


(b) fR\y=0 = 1,50 


Fig. 10 Plots of the effective potential W{z{y)) and relative probability P('w?) of the odd resonance KK modes on the brane 
generated by K-field with 6 = 1, d = —2, A = —0.01, and fji\y—Q = 1, 50. 


W(y) P 




Fig. 11 Plots of the effective potential W{z{y)) and relative probability P('m?) of the odd resonance KK modes on the brane 
generated by K-field with 6 = 8, d = —8, A = —0.01, and fji\y—Q = 1, 50. 


tential W{z{y)) of the /(i?)-brane on the whole space. 
For the cases with more scalar fields (such as the Bloch 
brane), one may draw the same conclusion, as long as 
all kinetic terms of the scalar fields are canonical. How¬ 
ever, for the brane generated by a K-field, the function 
fniy) could approach zero at some locations. There¬ 
fore, the effective potential W{z{y)) in the /(i?)-brane 
model may have plentiful structures and result in a se¬ 
ries of gravitational resonances, which is very different 
from the GR-brane model. 
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